This Page Is Inserted by IFW Operations 
and is not a part of the Official Record 



BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of 
the original documents submitted by the applicant. 

Defects in the images may include (but are not limited to): 

• BLACK BORDERS 

• TEXT CUT OFF AT TOP, BOTTOM OR SIDES 

• FADED TEXT 

• ILLEGIBLE TEXT 

• SKEWED/SLANTED IMAGES 

• COLORED PHOTOS 

• BLACK OR VERY BLACK AND WHITE DARK PHOTOS 

• GRAY SCALE DOCUMENTS 



IMAGES ARE BEST AVAILABLE COPY. 



As rescanning documents will not correct images, 
Please do not report the images to the 
Image Problem Mailbox. 




XP 000081157 <^ QfR /5g- A 

MAGNETIC RESONANCE IN MEDICINE 12 , 99-1 13 ( 1 989 ) 

8306 Magnetic Resonance in Medicine QC' R 33 f&&C2 
12(1989)0ctober, no.l, Duluth, MN, US ' 

Compensation for Effects of Linear Motion in MR Imaging 

Hope W. Korin,* Farhad Farzaneh,* Ronald C. Wright,* 
and Stephen J. Riederer 

Department of Diagnostic Radiology. Mayo Clink. Rochester. Minnesota 55905 

Received July 1 3. 1988: revised December 1 5. 1988 

Various compensation methods for different types of motion during MR image acqui- 
sition have been proposed. Presented here is a postprocessing scheme for eliminating 
artifacts due to linear, intra-slice motion of known velocity. The data for each phase 
encoding or "view" acquired from a moving object are shown to differ from those which 
would be measured from the stationary object by a phase factor which depends on the 
object's displacement from a reference point. This derivation is then used to propose a 
correction scheme for linear motion in which all phase encodings measured at different 
positions of the moving object are collapsed onto the same reference position. After sub- 
sequent reconstruction, the object appears perfectly "focused." By selection of different 
reference positions, the method permits positioning of the compensated object as desired 
within the field of view of the image. This property allows the method to be extended to 
create sequences of corrected images with smooth object motion between frames of the 
sequence. The basic correction scheme and its variations were tested experimentally in 
phantom studies with velocities as large as 8 cm/s. e iw Acmemie pres*. inc. 



INTRODUCTION 

Motion of a subject during MR image acquisition causes major artifacts in the 
reconstructed images. It is a challenge not only to compensate for or eliminate such 
artifacts but also to categorize and explain them in the first place. For standard two- 
dimensional Fourier transform ( 2DFT ) imaging (7. 2) various methods of compen- 
sation have been devised. These include reordering the acquisition of the phase en- 
codings to address the problems of respiratory motion (J) and gradient moment null- 
ing to ensure that moving spins are not dephased (4, 5). 

The problem we wish to address in this work is that of correction for linear, intra- 
slice motion which can occur throughout the duration of a scan. Before discussing 
this in detail, however, we wish to differentiate between two classes of inconsistencies 
in the acquired data, both of which are due to motion and both of which cause arti- 
facts in the reconstructed images. As in the distinction made by Pattany et at. (4), 
we call these 4t intra-view" and "intcr-view" effects. Here the term view refers to the 
process of the measurement of a single phase encoding. Intra-view inconsistency is 
caused by motion of an object in the interval between the initial RF pulse and the 
readout of the echo. A common example is the in-plane motion of blood with a 
range of velocities. Because of velocity dispersion, there will be a difference in the net 
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accumulation of phase bv spins moving ai different velocities, and upon reconstruc- 
tion some blood sienals will be reconstructed at only a fraction of their ideal value. 
Because of the manner in which the frequency-encoding gradients are applied, this 
phenomenon is eliminated in the second echo of a two-echo spin-echo sequence (6) 
if equal temporal spacing between echoes is chosen. Another example of intra-view 
inconsistency is the misplacement of objects in the image due to motion which occurs 
between the phase-encoding and the readout gradients ( ."). Such effects can be re- 
duced in severitv if the total time spanned by the gradient signals is minimized: i.e.. 
if the phase encodine. for example, is applied immediately before the signal readout 
(8) Phase problems arc compounded if the net phase accumulation dirters trom view- 
to view such as if pulsatile blood How occurs during a scan (9). However, the prob- 
lem of differential phase accumulation can he substantially eliminated with gradient 
moment nulling ( 4. 5). which ensures that spins moving with linear or higher order 
motion have zero net phase at the first echo. 

Inter-view effects arc taken to be inconsistencies due to motion between views. A 
paradigm is rcspiratorv motion which, if uncorrected, leads to ghosting in the phase- 
encoded direction ( 10). A reordering of the acquisition of the phase encodings using 
ROPE (i. / / ) or Exorcist techniques can suppress the ghosts, but this does not elimi- 
nate blurring. 

In this manuscript, we consider another inter-view effect, that of linear motion 
within the imace slice w ith velocity of known direction and magnitude. Because the 
object is movine. its position is different from view to view. Thus^ihis situation is 
distinct from the blood velocity case mentioned above in which differential phase is 
accumulated within a view. The situation addressed here is an inter- view inconsis- 
tency while gradient moment nulling addresses an intra-view inconsistency. 

An investigation into the mathematics orthis situation not only produces an expla- 
nation for the motion artifacts observed but also reveals a method ot compensation 
bv postproccssine of the acquired image data. This compensation scheme can gener- 
ate a corrected image of the object where the exact position to which it « recon- 
structed within the field of view can be chosen as desired. Moreover, as we will show, 
it becomes possible to generate a continuous sequence of focused images of the mov- 
ing object. ... 

Clinically this scheme would be useful in imaging any object which moves within 
the slice during a scan. In particular, it could be used in conjunction with a fast- 
scanning sequence such as MR fluoroscopy ( 12) to produce images of high enough 
aualitv to enable localization of structures in real time during patient positioning. In 
this case the patient table could be translated through the gantry with a velocity under 
operaior control. Further, the scheme presented here can potentially be "tended to 
the important physiological case in which the in-slice motion is nonlinear. We discuss 
how such a scheme could be implemented under Discussion. 

In the following sections, we present the theory for signal alteration due to linear 
motion a compensation scheme for it. and experimental results. The experiments 
presented consist of correction of single images and generation ol continuous image 
sequences of the moving object. The latter incorporates both correction for motion 
artifacts and selective placement of the object within the field of view. 
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THEORY 

The MR signal measured during a standard spin-warp image acquisition can be 
described by (1. 2) 

s„(t) - jj '"<a\ y)e-' ya » r '>e-' y0 *- u dxd\\ [1 ] 

where a* is the frequency-encoding (or readout) direction, r is the phase-encoding 
direction. (7 v and G n arc the respective gradients. m{ \. y) is the transverse magnetiza- 
tion of the object in spatial coordinates. 7 is the gyromagnetic ratio, and is the 
duration of the y gradient, n is an index over the repetition number of the acquisition. 
It runs from 0 to A * - 1 , where A* is the total number of phase encodings. 

If we consider an object moving in the plane of the selected slice, we can describe 
its velocity by two orthogonal components. i\ and t\. in the frequency and phase 
directions. Its position in the .vand r coordinates at the lime of the measurement of 
the signal can then be described as <.v 0 - v 0 - ty„). w here is the time of mea- 
surement of the nth phase encoding and .v 0 and r 0 describe an arbitrary reference 
position, and i\r„ are the displacements of the object from this original reference. 
We are ignoring the time delays due to the digitization during each signal measure- 
ment because they are very much smaller than those due to the repetition times be- 
tween acquisition of different signals. Translating this into mathematics and inserting 
it into Eq. [ 1 ] gives a formula for the signal y„( / ) measured for the case of motion: 

s:,(t) = jj mix - lv„, y-vj n )e-^^e- 1 ^' "d.xdy. [2] 

Now. using the substitution of variables 

« - a - iy„ [3a] 
i = y-v x t n [3b] 



we find 



s' n {t) - jj mU. ,,)c-'>^^'.».r-'> <; . << *^'- ,, ^/i? [4a] 
= jj m(l rjU'-'^v*^^ [4b] 

= e -Mjy*' >e -'y<i<w jj mi £, 17 )e' nG ^^e'' n0 ^ d^di) [4c] 

= e^Wve'^s^Snii), [4d] 

where the last result follows from Eq. [ I ] . 

We note that j*„(/), the signal measured from the moving object, is simply 5„(/). 
the signal from the stationary object measured for that phase encoding, multiplied 
by a phase factor which depends on the velocity of the object and the lime at which the 
signal was measured. In other words, the phase factor depends on the displacement of 
the object at the time of the measurement from some arbitran reference point. 
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Fig. I. Motion in the phase-encoded direction ( uncompensated images on the left, compensated on the 
right ): (A ) static phantom. ( B ) velocity of I cm/s. (C) 2 cnVs. <D)4.2cm/s. and ( E) 8.4 cm/s. The two 
white bars in ( E) indicate the distance traveled by the phantom during image acquisition. The motion was 
in the upward direction. 



Rather than work with Eq. [4d] in its present form with its explicit dependence on 
the values of the gradients, we note that the phase factors can be convened to field of 
view ( FOV ) dependence as follows. 




Fig. 1 — Continued 
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If the phase encodings (views) are applied in a standard sequential manner, then 
we have 



AC,. 



[5] 



Recall that .V is the total number of views per image and n is an index over these views, 
running from 0 to N- 1 . AC , is the incremental change in the encoding gradient from 
view to view. In general. AC, is related to the field of view ( FOV) by ( 2) 



AC,= 



2* 



7', (FOV) 



Therefore, the first phase factor of Eq. [4d] can be rewritten as 

e -i 2wi n - 1 .V - 1 1/ : )r, /„/ FOV 

Similarly, the gradient in the .v ( readout ) direction can be written as 

lit 



Cv = 



>A/(FOV) 



[6] 



[7] 



[8] 



where A/ is the interval between the samples taken of the MR signals and we have 
assumed the same FOV in .ras that in y. We next note that t for the readout direction 
is equal to zero at the center of the echo signal. 



-(,„-i^. 



[9] 





Fig. 2. Motion in the frequency-encoded direction (uncompensated image on the left, compensated on 
the right), velocity of I cm/s. The motion was in the upward direction. 
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Fig. 3. Comparison of two compensated images in which the phantom moved an equal distance during 
the image acquisition lime. The image on the left was acquired with a velocity of 1 cm/ s. 7K of 35 ms. 
that on the right with a velocity of 2 cm /s. 71* of 1 7. 5 ms. 



where m is the index of the readout point and \f is the total number of readout points, 
similar to n and ;V in the phase-encoding direction. We can then rewrite Eq. [4d] as 
*'„(/) = ^-^«"-i^-n/:) S vFov e -.2rt«HAi-n/2»v./Fov jB(/ j > [ l0 ] 

This is the equivalent of Eq. [ 4d ] , but with the FOV rather than gradient dependence 
within the phase factors. 

It is thus apparent from Eq. [10] that it is possible to recover the signal of the 
stationary object by multiplying the data from the moving object by the inverse of 
these phase factors, namely 

Sti ( t ) = e ^^^' i ^-^r.n\tjfov c ^r,wim- { \t-i)/Z)xjjroy s ' t{l y [\\] 

This is the basic correction scheme. For the cases of separate motion in either the 
frequency or phase-encoding direction, Eq. (1 1 ] simplifies because i\ and r are re- 
spect ivclv zero in these cases. 

In Eq. [1 1 ]. a,,, as previously discussed, represents the displacement of the object 
from its stationary position. This displacement will be different for every view ac- 
quired, as the object moves a distance t>TR between each two consecutive views. 
Multiplying each phase encoding by its specific phase factor in essence collapses all 
of them onto the position chosen as the stationary reference. It is thus possible to 
create images of the object from the same data set in which the object is centered 
at different locations within the field of view simply by choosing different reference 
positions for it. 

EXPERIMENTAL RESULTS 

Single Images 

The compensation technique of Eq. [1 1 ] was implemented experimentally on a 
GE Signa MR system. A limited flip angle pulse sequence ( 13) with the "GRASS" 
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Fig. 4. Illustration of object positioning. In (A), the reference view is -64: in( B). the reference view is 
+63. Both images were generated from the same image data. Note the difference in object position. 

variation ( 14) was used with TR/TE values of 17.5/8.9 ms and a field of view of 24 
cm. The object was moved by placing it on a table and translating the table through 
the gantry at an operator-controlled speed. A slice transecting the end of the phantom 
and parallel to the tabletop was imaged. By adjusting the pulse sequence, phase or 
frequency encoding was applied in the direction of the motion as desired. One hun- 
dred twenty-eight phase encodings were applied sequentially from most negative to 
most positive and then repeated. 

-64, -63 -1,0, +1 +63, -64', -63' 

the prime indicating the second application of a given encoding. Images were recon- 
structed from 128 consecutive views, and successive images were displaced by 32 
encodings. That is. Image 1 was formed from views -64 to +63. Image 2 from -48 
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Fig. 4— Continued 



to +63. and -64' to -49', etc. This is in essence the MR fluoroscopy reconstruction 
described in Ref.(/2). 

The first study is presented in Figs. 1 to 3. Figure 1 shows an image of the stationary 
object ( A ) . as well as uncompensated and compensated images of the object moving 
in the phase-encoded direction at speeds of 1 (B), 2 (C) t 4(D), and 8(E) cm/s. In 
all the images shown, the motion was in the upward direction. Figure 1 D shows a 
slightly different slice through the phantom than the other images. In Fig. 1 E, the two 
white lines at the top and near the bottom of the compensated image indicate the 
extremes of the distance traveled by the crossbar of the phantom during the acquisi- 
tion of the 128 views necessary to create the image. Figure 2 shows uncompensated 
and compensated images of the object moving in the frequency-encoded direction at 
a speed of I cm/s. Again the motion was upward. It is interesting to note the differ- 
ence in the appearance of the artifacts between motion in the frequency-encoded 
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Fic. 5. Schematic of creation of image sequences with MR fluorp' copy. ( A ) Object position with no 
compensation scheme implemented. ( B) continuous motion possibl *uh compensation. 

direction (Fig. 2, left) and phase-encoded direction (Figs. IB- IE. left). The slight 
curvature of the periphery of the image is due to spr .ial distortion caused by gradient 
nonlinearitv. However, this does not affect the fo.using of the image. Figure 3 is a 
comparison of two compensated images in whicl: the phantom moved an equal dis- 
tance during the image acquisition time. The ir .age on the left was acquired with a 
velocity of 1 cm/s and a TR of 35 ms, and that on the right with a velocity of 2 cm/s 
and a TR of 17.5 ms. That is. in each case tite incremental displacement between 
successive phase encodings was uTR or 0.35 .nm. As expected, the quality of the two 
images is very similar. The pattern of the ar.ifacts in the uncompensated images ( not 
shown) is also similar. ... 

In all the images of Figs. 1 to 3. the middle of the time interval of data acquisition 
for the image was chosen as the reference position for the object. As can be seen, each 
of the compensated images is of quality comparable to that of the stationary image 
( Fig 1 A ) 

The second studv involved investigating the modification of a single data set to 
generate images with various object positions. A series of images was created, all from 
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Fig. 7. Illustration of velocity tuning. The actual velocity was 8.4 cm/s. shown on the rieht. The image 
on the left shows the effect of implementing the compensation scheme on these data assuming the incorrect 
velocity of 8.2 cm/s. 

the same data set of 1 28 views ( numbered -64 to +63 ). such that in each image the 
measurement time of a different phase encoding was selected as the reference point. 
This is equivalent to choosing the time at which the reference phase encoding was 
measured to be the zero of the time scale. Figure 4 shows two of these compensated 
images generated from the same data set with phantom motion of 1 cm/s in the 
phase-encoding direction. In Fig. 4A, the reference view is -64. the first encoding of 
128 acquired. That is, although each view was measured with a different object posi- 
tion, each was multiplied by the appropriate phase factor to shift its object position 
back to that measured by the -64 view. In Fig. 4B. the reference encoding is +63, 
the last one measured. Here the views arc collapsed forward onto it. The difference 
in position of the object between the two images is easily seen. 



Image Sequences 

This technique is of great interest for application to MR fluoroscopy imaging ( 12). 
Fluoroscopy is a fast-scanning technique that seeks to produce images of dynamic 
systems in real time, such as in patient positioning before a scan. In this case, a se- 
quence of corrected images of the moving object would be desirable to produce a 
video sequence of the motion. The two main undesirable effects observed in uncom- 
pensated videos of moving objects are the motion artifacts themselves and the large 
"hops" or "jumps" of object position seen instead of continuous motion between 
video frames ( 12). In some frames, a duplication of the object is even observed. An 
explanation for the large hopping movements and the object duplication is shown in 
Fig. 5 A. This figure is a plot versus time of the application of phase encodings. One 
can also imagine for purposes of discussion that the object is moving to the right. 
Two repetitions of the application of 123 views are shown, as discussed previously. 
Images are formed by taking a set of 1 28 contiguous phase encodings. Here the offset 
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between successive images is chosen to be 16 views for illustrative purposes. The low 
spatial frequencies represent the encoding of the gross position of the object. Higher 
frequencies represent the edges. Let us assume initially that no motion compensation 
is done. Then Images 1 through 4 will show the gross object position to be the same 
since they all use the same low-frequency measurements in their reconstruction. It is 
the edges of the object which will be updated in going from Image I to 4 and in fact, 
the edges of the object appear to move ahead of the object. Image 5 will show the 
effect of duplication of the object because it incorporates low spatial frequency mea- 
surements from two distinct times and thus two distinct object positions. Image 6 
will show the second gross position of the object only and so forth. Because the gross 
position of the object will change only when a new set of low frequencies is measured, 
it appears to be stationary for a number of images (4 or 5 in this illustration) and 
then will hop a considerable distance on the next image. 

It is possible to eliminate these effects by applying the features of the correction 
scheme discussed previously under single image compensation. This is illustrated in 
Fig. 5B. By implementing the correction method, not only are the motion artifacts 
eliminated, but also the hopping motion is eradicated by a judicious choice of the 
reference position of the object in each image. If we choose the position of the object 
at the center of the data acquisition time to be the reference position of the object 
(i.e.. view 0 in Image 1. view 16 in Image 2. and view 32 in Image 3) we can generate 
a sequence of images of the moving object in which the distance traveled by the object 
between successive frames is always the same: thus a clear, continuously moving ob- 
ject is seen. 

Figures 6A and 6B show stills of such an image sequence. This sequence w as recon- 
structed with a displacement between images of 32 views instead of the 16 show n in 
the drawings of Fig. 5. but the principle is the same. The center portions of four 
uncompensated images from the sequence are shown in (A), the corresponding com- 
pensated ones in ( B ). It should be stressed that the same data set was used for each 
sequence. Image 1 (panel ! in (B). numbering the panels from left to right) was 
compensated to view 0. Image 2 (panel 2) to view 32, Image 3 (panel 3) to view -64. 
and Image 4 (panel 4) to view -32. Note the apparent continuous motion of focused 
images in (B) as opposed to the hopping motion in the uncompensated images of 
( A ) and the elimination of the duplication of the object in panel 3. 

DISCUSSION 

In the work presented here, we have described how the signals acquired during a 
spin-warp pulse sequence are modified if an object moving at constant, linear velocity 
within the slice is imaged. It was shown that such signals differ from the signal that 
would be acquired from a stationary object by a phase factor, different for each view 
and dependent on the distance the object has traveled from some reference position. 
It is thus possible to generate an uncorrupted image of the object by multiplying the 
data acquired by the inverse of these phase factors. Moreover, as there is no prede- 
fined reference position, it is possible to choose it as desired, and thus position the 
object arbitrarily within the field of view during the postprocessing scheme. 

Because the experiments described here were performed with a moving phantom, 
one might question whether there were any intra-view inconsistencies in the acquired 
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data Because the entire object was moving at the same velocity, the accumulated 
phase of all spins within a single view was identical. Additionally, because the velocity 
was nominally constant throughout the entire scan, this phase accumulation was the 
same for every view. Hence, there were no inconsistencies due to intra-view motion. 
The effects were solely due. as desired, to different vicw-to-view displacements of the 
object from a reference position. 

The basic correction scheme was verified experimentally with velocities as high as 
8 cm / s. At this speed, the object movement during acquisition of all the views neces- 
sarv to reconstruct one image was about 75<?c of the field of view. The modified tech- 
nique of object placement within the field of view was also investigated. Both methods 
were then used to generate sequences of images in which the motion from one video 
frame to the next was not only free of artifacts, but also displayed smooth and contin- 
uous motion. . . 

A potential extension to this method is velocity tuning, in which the operator could 
retrospectively determine the velocity at which an object was traveling by adjusting 
the phase factors until the images cleared. An example of implementing the correc- 
tion scheme with a slightly incorrect velocity is shown in Fig. 7. The image 
on the left was reconstructed assuming a velocity of 8.2 cm/s. instead of the correct 
8 4 cm/s shown on the right. Note especially the focusing of the lines m the vertical 
crossbar of the phantom. Clearly, if one attempts to image a collection of objects 
which are moving at different rates, this method alone would not permit, focusing of 
all of them at once. Also, parts of the object entering the slice for the first time may 

not be corrected properly. c ^;,i 
In the derivation and experimental results desenbed here, we assumed the special 
case of linear motion with known velocity. Other than the MR fluoroscopy /moving 
table example mentioned earlier, one may question how this can relate to clinical 
imaging. In fact, the derivation and correction can easily be modified to allow Tor a 
more general case. Suppose in the nth view that the displacement of the object of 
interest from a reference position is some arbitrary .v„ rather than the r,/„ assumed in 
Eq [21 In this case, it is straightforward to show that the same substitution (x n for 
vj m ) is made in the correction term of Eq. [4d]. Having derived this, one next ques- 
tions how the set of displacements {.x„} can be measured. In fact this can be readily 
accomplished using a dual-echo pulse sequence in which one echo encodes the de- 
sired image information and the second echo (a "monitor" or -navigator echo at 
the zeroth phase encoding) is used to monitor the view-to-view displacements of a 
potentially moving structure ( 15 ) . This concept of using a navigator echo in conjunc- 
tion with the data correction scheme presented in this work has been proposed by 
Ehman el ai ( 16). It is expected that such a technique can more completely correct 
clinical images for motion blurring than present methods. 
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